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ABSTRACT 

We present broad-band radio observations of the afterglow of GRB 00030 1C, spanning from 1 .4 to 350 GHz for 
the period of 3 to 83 days after the burst. This radio data, in addition to measurements at the optical bands, suggest 
that the afterglow arises from a collimated outflow, i.e. a jet. To test this hypothesis in a self-consistent manner, we 
employ a global fit and find that a model of a jet, expanding into a constant density medium (ISM+jet), provides 
the best fit to the data. A model of the burst occurring in a wind-shaped circumburst medium (wind-only model) 
can be ruled out, and a wind+jet model provides a much poorer fit of the optical/IR data than the ISM+jet model. 
In addition, we present the first clear indication that the reported fluctuations in the optical/IR are achromatic with 
similar amplitudes in all bands, and possibly extend into the radio regime. Using the parameters derived from the 
global fit, in particular a jet break time, tj e , w 7.5 days, we infer a jet opening angle of 9q ~ 0.2, and consequently 
the estimate of the emitted energy in the GRB itself is reduced by a factor ^50 relative to the isotropic value, 
giving £»l,lx 10 51 ergs. 

Subject headings: gamma rays:bursts - radio continuum:general - cosmology:observations 



1. INTRODUCTION 

GRB 00030 1C is the latest afterglow to exhibit a break in its 
optical/IR light curves. An achromatic steepening of the light 
curves ha s been interpreted in previous events (e.g. Kulkarni et 
al. 1999 a; Harrison et al. 1999[) as the signature"o f a jet-like 
outflow flRhoads 1999| ; [Sari Piran & Halpern 1999| ), produced 
when relativistic beaming no longer "hides" the non-spherical 
surface, and when the ejecta undergo rapid lateral expansion. 
The question of whether the relativistic outflows from gamma- 
ray bursts (GRBs) emerge isotropically or are collimated in jets 
is an important one. The answer has an impact both on esti- 
mates of the GRB event rate and the total emitted energy — 
issues that have a direct bearing on GRB progenitor models. 

An attempt by Rhoads and Fruchter (2000) to model this 
break using only the early time (Af < 14 days) optical/IR data 
has led to a jet interpretation of the afterglow evolution, but 
with certain peculiar aspects, such as a different jet break time 
at R band than at K' band. However, subsequent papers by 
Masetti et al. (2000), and Sagar et al. (2000), with larger optical 
data sets, pointed out that there are large flux density variations 
(~ 30%) on timescales as short as a few hours, superposed on 
the overall steepening of the optical/IR light curves. While the 
origin of these peculiar fluctuations remains unknown, it is clear 
that they complicate the fitting of the optical/IR data, rendering 
some of the Rhoads and Fruchter results questionable. 

In this paper we take a different approach. We begin by 



presenting radio measurements of this burst from 1 .4 GHz to 
350 GHz, spanning a time range from 3 to 83 days after the 
burst. These radio measurements, together with the published 
optical/IR data, present a much more comprehensive data set, 
which is less susceptible to the effects of the short-timescale op- 
tical fluctuations. We then use the entire data set to fit a global, 
self-consistent jet model, and derive certain parameters of the 
GRB from this model. Finally, we explore the possibility of 
a wind, and wind+jet global fit to the data, and compare our 
results with the conclusions drawn in the previous papers. 

2. OBSERVATIONS 

Radio observations were made from 1 .43 GHz to 350 GHz, at 
a number of facilities, including the James Clark Maxwell Tele- 
scope (JCMT 8 ), the Institut fur RadioAstronomie im Millime- 
terbereich (IRAM 9 ), the Owens Valley Radio Observatory In- 
terferometer (OVRO), the Ryle Telescope and the (Very Large 
Array (VLA 10 ). A log of these observations and the flux den- 
sity measurements are summarized in Table [l] With the excep- 
tion of IRAM, we have detailed our observing and calibration 
methodology in Kulkarni et al. (1999) and Frail et al. (2000). 

Observations at IRAM were made using the Max-Planck 
Millimeter Bolometer (MAMBO; Kreysa et al. 1999) at the 
IRAM 30-m telescope on Pico Veleta, Spain. Observations 
were made in standard on-off mode. Gain calibration was per- 
formed using observations of Mars, Uranus, and Ceres. We es- 
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timate the calibration to be accurate to 15%. Using the MOPSI 
software package (Zylka 1998), the temporally correlated vari- 
ation of the sky signal (sky-noise) was subtracted from all 
bolometer signals. The source was observed on March 4, 5, 
and 9 under very stable atmospheric conditions, and on March 
6 with high atmospheric opacity. From March 24 to 26, the 
source was briefly re-observed three times for a total on+off 
integration time of 2000 sec, but no signal was detected. 

3. THE DATA 

In Figure [j] we present broad-band spectra from March 5.66 
UT (Af ps 4.25 days) and March 13.58 UT (At « 12.17 days). 
Radio light curves at 4.86, 8.46 22.5, and 250 GHz from Ta- 
ble U are presented in Figure |J while optical/IR curves are 
shown in Figure ^| 

The quoted uncertainties in the flux densities given in Ta- 
ble |l] report only measurement error and do not contain an es- 
timate of the effects of interstellar scatterin g, which is know n 
to be significant for radio afterglows (e.g. Frail et al. 1999 ). 
We can get some guidance on the expected magnitude of the 
IS S -induced modulation of our flux density measurements (in 
time and frequency) using the models developed by Taylor and 
Cordes (1993), Walker (1998), and Goodman (1997). 

From the Galactic coordinates of GRB 00030 1C 
(/,6)=(48.7°,44.3°), we find, using the Taylor and Cordes 
model, that the scattering measure is 5M_3 5 ps 0.7. The dis- 
tance to the scattering screen, d scr , is one half the distance 
through the ionized gas layer, d scr = (hz/2)(sinb)~ 1 ps 0.72 kpc, 
using hz ~ 1 kpc. From Walker's analysis, the transition fre- 
quency between weak and strong scintillation is then given by 
u = 5.9SM%"d 5 J r 11 « 4.7 GHz. Goodman (1997) uses the 
same scalings, but with a different normalization for the tran- 
sition frequency, giving a larger value, vq ps 8.3 GHz. In this 
section we follow Walker's analysis, and note that the numbers 
from Goodman will give different results. 

For frequencies larger than the transition frequency the mod- 
ulation index (i.e. the r.m.s. fractional flux variation) is 
m v = (uq/i/) 11 ! 12 , and the modulation timescale in hours is 
?„ ps 6.7(d scr / V) 1 / 2 . From this analysis we find that the mod- 
ulation index is of order 0.4 at 8.46 GHz, 0.2 at 15 GHz, 0.10 
at 22.5 GHz, and is negligible at higher frequencies. The mod- 
ulation timescales are of order 2.0 hours at 8.46 GHz, 1.5 hours 
at 15 GHz, and 1.2 hours at 22.5 GHz. It is important to note 
that factor 2 uncertainties in the scattering measure allow the 
modulation index to vary by ~ 50%. 

At these frequencies the expansion of the fireball will be- 
gin to "quench" the ISS when the angular size of the fire- 
ball exceeds the angular size of the first Fresnel zone, Of = 
8(af JC/ .^GHz)~ 1//2 Ai as - To describe the evolution of the source size 
with tim e, we have used an expanding jet model (see F rail et 
al. 1999), with the factor (E^/ni) 1 ^ assumed to be of order 
unity, which gives 9 S ps 3.1(Af^/15) 1 ' 2 /ias. Once the source 
size exceeds the Fresnel size (after approximately two weeks at 
8.46 GHz), the modulation index has to be corrected by a factor 
(Au/IST 1 ' 11 . 

The measurements at 4.86 GHz occur near the transition 
frequency and we therefore expect 1114 gg to be large ~ 0.65 - 
1.0. At 1.43 GHz the observations were made in the strong 
regime of ISS where we expect both refractive and diffrac- 
tive scintillation. Point source refractive scintillation at 1.43 
GHz has a modulation index mi.43 r = (v / vq) 11 / 30 ps 0.5, with a 
timescale of ?i.43, r ps 2(vo/v) n l 5 ps 1 day. The refractive ISS is 



"quenched" when the angular size of the source is larger than 
9 r = 6F()(va/v) n l 5 , where 6po is the angular size of the first 
Fresnel zone at v$ = 4.7 GHz. As with weak scattering, the 
modulation index has to be corrected by a factor (Afj/15) -7 / 12 
after this point. The diffractive scintillation has a modulation 
index nix ^j = 1, and a timescale, fi.43^ ps 2(v / Vo) 6//5 ps 0.5 
hrs <C fi.43, r . The source can no longer be approximated by 
a point source when its angular size exceeds 9d = Ofo^/vo) 6 ^, 
and correspondingly the modulation index has to be corrected 
byafactor(Af rf /15)- 1 / 2 . 

The redshift of GRB 00030 1C was measured using the Hub- 
ble Space Telescope to be 1.95 ±0.1 by Smette et al. (2000) 
and was later refined by Castro et al. (2000) using the Keck II 
10-m Telescope to a value of 2.0335±0.0003. The combined 
fluence measured by the GRB detector on board the Ulysses 
satellite, and the X-ray/gamma-ray Spectrometer (XGRS) on 
board the Near Earth Asteroid Rendezvous (NEAR) satellites, 
in the 25 - 1 00 keV and > 1 00 ke V bands was -4 x 1 0" 6 erg/cm 2 . 
Using the cosmological parameters f2o=0.3, Ao=0.7 and //q=65 
km/sec/Mpc, we find that the isotropic 7-ray energy release 
from the GRB was E iso ps 5.4 x 10 52 ergs. 

4. A SELF-CONSISTENT JET INTERPRETATION 

According to the standard, spherical GRB model, the optical 
light curves should obey a simple power-law decay, F v oc t~ a , 
with a changing at most by 1/4 as the electrons age and cool 
(Sari, Piran, & Narayan 1998). From Figure [| it is evident that 
the optical light curves steepen substantially (Aq > 1/4) be- 
tween days 7 and 8, which indicates that this burst cannot be 
described within this standard model of an expanding spherical 
blast wave. This break can be attributed to a jet-like or colli- 
mated ejecta (Rhoads 1999; Sari, Piran, & Halpern 1999). 

The jet model of GRBs predicts the time evolution of flux 
from the afterglow, and of the parameters v a oc f -1 / 5 , u m oc t~ 2 , 
and F vmax oc t . This model holds for t > tj et , where tj et is 
defined by the condition j(tj et ) ~ 8q 1 . Prior to tj et the time 
evolution of the afterglow is described by a spherically ex- 
panding blast wave, with the scalings v a oc const., v m oc f~ 3//2 , 
and F umax oc const. In this paper we designate this model as 
ISM+jet. Throughout the analysis we assume that the cooling 
frequency, i/ c , lies above the optical band for the entire time 
period under discussion in this paper. 

At any point in time the spectrum is roughly given by the bro- 
ken power law F v oc v 1 for v < u a , F v oc i/ 1 / 3 for v a < v < u m , 
and F v oc v~ (p ~ 1 ^ 2 for v > v m , where p is the electron power 
law index. To globally fit the entire radio and optical/IR data 
set we employed the smoothed form of the broken power law 
synchrotron spectrum, calculated by Granot, Piran, and Sari 
(1999a,b). With this approach we treat tj et , p, and the values 
of v a , v m , and F vmax at t = tj et as free parameters. This method 
forces tj e t to have the same value at all frequencies. We find 
the following values for the burst parameters: tj e , = 7.5 ±0.5 
days, p = 2.70 ± 0.04, v a (t = t jet ) = 8.3 ±1.8 GHz, v m (t = t jet ) = 
(3.0 ±0.3) x 10 11 Hz, F v . mx (t = t jet ) = 2.7 ±0.2 mJy, where 
the errors are derived from the diagonal elements of the cor- 
relation matrix. We note that there is substantial covariance 
between some of the parameters and therefore these error esti- 
mates should be treated with caution. From our fit the asymp- 
totic temporal decay slopes of the optical light curves are p we 
find ai = -3(/?- l)/4 = -1.28 for t < t je „ and a 2 = -p = -2.70 
for t > tjet ■ The fits are shown in figures [IJ |^, and ||. 

The total value of x 2 for the global fit is poor. We obtain 
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X 2 = 450 for 94 degrees of freedom. The bulk of this value, 340, 
comes from the 61 optical data points, and is the result of the 
fluctuations, which are not accounted for by our model. The ra- 
dio data contribute a value of 1 10 to % 2 for 38 data points. This 
is probably the result of scintillation. If we increase the errors 
to accommodate for the expected level of scintillation (see §3) 
we obtain a good fit with x\adw = 39/33 degrees of freedom. 

From figure |l] it is clear that the global fit accurately de- 
scribes the broad-band spectra from days 4.26 (< tj et ) and 12.17 
(> tj e t), with a single value of p = 2.70, which rules out the pos- 
sibility that the steepening of the light curves at t = tj et is the 
result of a time- varying p. 

Trying to model the data using the approach outlined above, 
but for a wind-shaped circumburst medium results in a poor de- 
scription of the data, because the wind model does not exhibit 
a break, although one is clearly seen in the optical data. As a 
result, the model fit is too low at early times, and too high at 
late times relative to the data (see insert in figure ||). The value 
of x 2 f° r the wind model relative to the ISM+jet model de- 
scribed above is xtirui/XisM+jet ~ 3. Therefore, a wind-shaped 
model can be ruled out as a description of the afterglow of 
GRB000301C. 

A jet evolution combined with a wind-shaped circumburst 
medium provides a more reasonable fit than a wind only model. 
The wind evolution of the fireball will only be manifested for 
t < tj e t since once ~f(tj et ) « 9q 1 the jet will expand sideways 
and appear to observers as if it were expanding into a con- 
stant density medium (Chevalier & Li 1999b; Livio & Wax- 
man 1999). The resulting parameters from such a fit differ con- 
siderably from the parameters for the ISM+jet model quoted 
above, and the relative value of x 2 between the two models is 
xlmd+jet/XisM+jet ~ 2 - This model suffers from a serious draw- 
back in its description of the optical/IR light curves. Because 
the predicted decay of these light curves prior to f /e , is steeper 
than in the ISM+jet model, the model fit, from 2 days after the 
burst up to the break time, is too low relative to the data (see 
insert in figure [|). 

The approach outlined above, which uses one power law tem- 
poral evolution of v a , v m , and F vmax before tjet, an d a different 
power law evolution after tj et , creates a sharp break at tj et as 
seen in Figures ||, and || In contrast, a smooth analytical form 
for the jet break, was used by several other groups (Harrison 
et al. 1999; Stanek et al. 1999; Israel et al. 1999; Kuulkers 
et al. 2000) to describe the afterglow of GRB 990510. Beuer- 
mann et al. (1999), suggested that since there is no detailed 
theory for the jet transition, the shape of the break (i.e. how 
smooth or sharp it is) should be kept as a free parameter. In 
a recent paper, Rhoads and Fruchter (2000) used the same ap- 
proach, a free shape-parameter, and got a smooth light curve 
as their best fit. However, they have not forced the relation be- 
tween the asymptotic temporal decay slopes before and after 
tj e t, (c*i and a2, respectively) and they allowed separate slopes 
and break times for the R and K' bands. Using the Beuermann 
et al. formula, forcing the asymptotic relations between ot\ and 
OL2 to fit the theory, and using a single tj et for both bands, we 
find that the best fit is a sharp break. This may be the result 
of the unexplained fluctuations that appear in the optical bands. 
In the radio regime, there is no data around tj et and therefore a 
smooth connection would not make a difference. This justifies 
the use of a sharp break in our fit. 

The global fitting approach has several advantages over fit- 
ting each component of the data set independently. For exam- 
ple, the K' data is only available up to day 7.18 (^ tj et ) after the 



burst. Therefore, by fitting it independently of the R band and 
of the radio data we cannot find tj et , if it is indeed after 7 days. 
Moreover, since, as Masetti et al., and Sagar et al. claim, and as 
we can see from Figure |J there is an additional process which 
superposes achromatic fluctuations with an overall rise and de- 
cline centered on day 3, on top of the smoothly decaying optical 
emission (see insert in figure [|), then fitting the K' data inde- 
pendently will confuse this behavior with the jet break. This 
explains the result of Rhoads and Fruchter of tj e t,K' ~ 3 days. 
It is worth noting that fitting the available R band data from 
before day 8 by itself, gives a value of tj etl R ~ 3.5 days ~ tjetjc'. 

Simultaneous fitting of the entire data set makes it possi- 
ble to study the overall behavior of the fireball regardless of 
any additional sources of fluctuations, because the large range 
in frequency and time of the data reduces the influence of 
such fluctuations. Remarkably, using this global fit with only 
the radio data, ignoring the optical observations, we obtain 
tjet,Radio ~ 7.7days ~ tjet- Thus, the radio data serves to sup- 
port the jet model, and provides an additional estimate for the 
jet break time, independent of the somewhat ambiguous optical 
data. 

From the global fit we find the first self-consistent indica- 
tion that the short-timescale optical fluctuations are achromatic, 
even in the K' band (see insert in figure [}]). By simply divid- 
ing the B, R, V, I, and K' data by the values from the global 
fit we find that the fluctuations happen simultaneously and with 
similar amplitudes in all bands. Moreover, the overall structure 
of the fluctuations is a sharp rise and decline centered on day 
4, and with an overall width of 3.5 days, which gives St/t ~ 1, 
where St is the width of the bump. The optical/IR data starts 
at day 1.5 lower by 25-50% than the model fit, then rises to a 
peak level of 50-75% relative to the model at day 4, and drops 
to the predicted level at about day 5, at which point it follows 
the predicted decline of the ISM+jet model. 

It is interesting to note that the 250 GHz data, which is not 
affected by ISS-induced fluctuations, also shows a peak ampli- 
tude approximately 70% higer than the model fit around day 
4 (see insert in figure ||). At the lower radio frequencies there 
are not enough data points to discern a similar behavior. More- 
over, at these frequencies it would have been difficult to disen- 
tangle such fluctuations from ISS-induced fluctuations in any 
case. The large range in frequency of this achromatic fluctua- 
tion, coupled with the similar level of absolute deviation from 
the model fit suggests that it is the result of a real physical pro- 
cess. 

It is possible to explain this fluctuation as a result of a non- 
uniform ambient density. The value of v m is independent of 

l li 

the ambient medium density, and since F vmax oc «/ , we ex- 
pect the flux at frequencies larger than v a to vary achromati- 

1/2 

cally, and with the same amplitude, F v oc n/ . For frequen- 
cies lower than v a we have to take into account the density 

3/5 

dependence v a oc n/ so that the flux will vary according to 

F v oc Fv Jnax h>~ 2 oc nj 7//1 °. This means that for frequencies lower 
than ^11 GHz we actually expect the flux to fluctuate down- 
ward at the same time that it fluctuates upward at higher fre- 
quencies. In practice, we don not have enough data around this 
time to confirm this behavior. In order to match the observed 
amplitude of the fluctuation ~ 80%, the ambient density has to 
vary by about a factor of 3. 

Using the value of t j er from our global fit, we can calculate 
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the jet opening angle, 9q, from the equation: 

9 w 0.05(f > //^) 3/8 d +zr 3/8 (ni/£ 5 2) 1/8 , (D 

(Sari et al. 1999; Livio & Waxman 1999) where £52 is the 
isotropic energy release, which can be roughly estimated from 
the observed fiuence (using the equations from Rhoads (1999) 
results in a smaller opening angle). From this equation, we 

1 /8 

calculate a value of 6*o ~ 0.2«j radians. This means that the 
actual energy release from GRB 00030 1C is reduced by a fac- 
tor of 50 relative to the isotropic value, Ei SO m 5.4 x 10 52 ergs, 
which gives £=1.1 x 10 51 nj ergs. 

5. CONCLUSION 

The afterglow emission from GRB 00030 1C can be well de- 
scribed in the framework of the jet model of GRBs. Global 
fitting of the radio and optical data, allows us to calculate the 
values of p, tj et , and the time evolution of v a , v m , and F Vimax 
in a self-consistent manner. Within this approach the proposed 
discrepancy between the behaviors of the R band and K' band 
light curves, suggested by Rhoads and Fruchter, is explained as 
the result of the lack of data for t > 7.18 days (^ tj et ) at K', 
and the existence of achromatic substructure from fluctuations 



in the radio and optical/IR regimes. The value for the break 
time from the global, self-consistent approach we have used is 
tj et = 7.5 days at all frequencies. 

The long-lived radio emission from the burst, spanning a 
large range in frequency and time, plays a significant role in 
our ability to extract the time evolution of v a , v m , and F Uttnax 
from the data. In the case of this GRB in particular, the large 
range in frequency and time is crucial, since it serves to reduce 
the effects of unexplained deviations from the simple theory, 
such as the short-timescale fluctuations in the optical bands, on 
the overall evolution of the fireball. 

We end with some words of caution. In our analysis we as- 
sumed that v c lies above the optical band throughout the evo- 
lution of the fireball, and successfully got a reasonable fit to 
the data. However, it is possible that another set of parameters, 
with v c below the optical band, can fit the data equally well. 
Preliminary work in this direction indicates that the gross fea- 
tures of the fireball evolution (e.g. a break time of ~7 days) 
remain unaltered. 

Research at the Owens Valley Radio Observatory is sup- 
ported by the National Science Foundation through NSF grant 
number AST 96-13717. 
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Table 1 

Radio Observations of GRB 00030 1C 



Epoch 


A t 


lelescope 






fl TT"\ 


(days) 








2000 March 4.29 


2.88 


IRAM 


250 


2100±300 


2000 March 4.75 


3.34 


JCMT 


350 


3736±3700 


2000 March 4.98 


3.57 


Ryle 


15.0 


660±160 


2000 March 5.41 


4.00 


IRAM 


250 


2300±400 


2000 March 5.53 


4.12 


JCMT 


350 


2660±1480 


2000 March 5.57 


4.16 


OVRO 


100 


2850±950 


2000 March 5.67 


4.26 


VLA 


1.43 


11 ±79 


2000 March 5.67 


4.26 


VLA 


4.86 


240±53 


2000 March 5.67 


4.26 


VLA 


8.46 


316±41 


2000 March 5.67 


4.26 


VLA 


22.5 


884±216 


2000 March 6.29 


4.88 


IRAM 


250 


2000±500 


2000 March 6.39 


4.98 


VLA 


8.46 


289±34 


2000 March 6.50 


5.09 


JCMT 


350 


1483±1043 


2000 March 6.57 


5.16 


OVRO 


100 


-99±1500 


2000 March 9.25 


7.84 


IRAM 


250 


400±600 


2000 March 10.21 


8.80 


Ryle 


15.0 


480±300 


2000 March 13.58 


12.17 


VLA 


8.46 


483±26 


2000 March 13.58 


12.17 


VLA 


22.5 


748±132 


2000 March 15.58 


14.17 


VLA 


8.46 


312±62 


2000 March 17.61 


16.20 


VLA 


8.46 


380±29 


2000 March 21.52 


20.12 


VLA 


8.46 


324±36 


2000 March 23.55 


22.14 


VLA 


8.46 


338±69 


2000 March 24.29 


22.88 


IRAM 


250 


-300±500 


2000 March 27.55 


26.14 


VLA 


8.46 


281±34 


2000 March 31.53 


30.12 


VLA 


8.46 


281±25 


2000 April 4.59 


34.18 


VLA 


8.46 


325±27 


2000 April 10.36 


39.95 


VLA 


8.46 


227±33 


2000 April 12.47 


42.06 


VLA 


4.86 


210±43 


2000 April 12.47 


42.06 


VLA 


8.46 


91±38 


2000 April 15.43 


45.02 


VLA 


8.46 


233±37 


2000 April 18.47 


48.06 


VLA 


4.86 


226±51 


2000 April 18.47 


48.06 


VLA 


8.46 


145±36 


2000 May 4.49 


64.13 


VLA 


4.86 


136±45 


2000 May 4.49 


64.13 


VLA 


8.46 


150±20 


2000 May 7.50 


67.09 


VLA 


4.86 


85±33 


2000 May 7.50 


67.09 


VLA 


8.46 


144±31 


2000 May 22.45 


82.04 


VLA 


8.46 


105±25 


2000 May 23.45 


83.04 


VLA 


8.46 


114±24 



NOTE. — The columns are (left to right), (1) UT date of the start of each observation, (2) time elapsed since the 7-ray burst, (3) telescope name, (4) observing 
frequency, and (5) peak flux density at the best fit position of the radio transient, with the error given as the root mean square noise on the image. The JCMT 
observations did not detect the source at each epoch individually, but by averaging the 3.875 hours of integration over the three epochs, we obtain a 2.5ct detection of 
1.70 ± 0.71 mJy. 
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FIG. 1 .— The radio to optical spectral flux distribution of GRB 000301C on 2000 March 5.66 UT (At ss 4.26 days after the burst), and 2000 March 13.58 UT 
(Af ?3 12.17 days). The solid and dashed lines are the global fits based on the smoothed synchrotron emission spectrum of Granot, Piran, and Sari (1999a,b). The 



radio mpasiirpmpnts arp from Tahlp |l| Thp nptipal/TR Hata arp from Rtinark and Fmchter (2000), Sagar et al. (2000), and Maytti pt al (2000), rnnvprtpH tn TangVy 

flux units ( Bessell & Brett 1988] , |Fiikugita Shimasaku & Ichikawa 199' ), and corrected for Galactic foreground extinction (Schlegel, Finkbeiner & Davis 1998), 
giving E(B- V) = 0.U53. All data were taken within 0.5 days of the fiducial dates, and the circles are the corrections to the data to the fiducial times, = 4.26 days, 
and A; = 12.17 days. The squares in the optical bands are weighted averages of multiple measurements within 1 day of Af = 4.26 days (see insert). The inverted 
triangle is an upner limit at 1.43 GHz from day 4.26. The data points at 100, 250, and 350 GHz are weighted averages of the individual measurements from around 
day 4 (see Table hi). Note that the data and fit from Af = 12.17 days were divided by a factor of 10 to avoid overlap with the Af = 4.26 curve. 
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FIG. 2. — Radio light curves for 4.86, 8.46, 22.5, and 250 GHz. The data were obtained at the VLA and IRAM (see §2). The solid lines are the global fit, based 
on the time dependences of the parameters u a , u,„, and F Vymax , and the smoothed, broken power law synchrotron spectrum calculated by Granot, Piran and Sari (see 
§4). The sharp break at time t = tj et ~ 7.5 days corresponds to the transition from spherical to jet geometry. The dashed curve shows the prediction for a spherical 
evolution of the afterglow. The dotted lines indicate the maximum and minimum range of flux expected from ISS (see §3). Note that the data and fit for 4.86 GHz 
were divided by a factor of 10, the data and fit for 22.5 GHz were multiplied by a factor of 10, and the data and fit for 250 GHz were multiplied by a factor of 100 
to avoid overlap between the four curves. 
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FIG. 3.— Optical light curves for B, R, and K' bands of GRB 000301C. The data are from Rhoadiand Fruchter (2000), Sagar et al. (2000), Masetti et al. (2000), 
and Diercks et al. (2000; filled circle), and contain a correction for Galactic extinction (see Figure hj). Following Masseti et al. (2000), we added a 5% systematic 
uncertainty in quadrature to all optical measurements to account for discrepancies between the different telescopes and instruments used in the observations. The 
solid lines are the global fit, based on the smoothed, broken power law synchrotron spectrum calculated by Granot, Piran and Sari (see §4). In the top right insert are 
plotted the data points divided by the respective model fit for all bands (circles, squares, crosses, triangles, pluses, and inverted triangles indicate R, B, K', V, I, and 
250 GHz bands, respectively). It is evident that the short-timescale fluctuations are achromatic and with a comparable amplitude in all bands, spanning from optical 
to radio. The insert on the bottom left portion of the figure shows the global fits based on the wind-only (dotted line), and wind+jet (dashed line) models overlaid 
on the R band data. It is clear that the steeper decline predicted for a fireball expanding into a wind-shaped circumburst medium results in a much poorer fit relative 
to the ISM+jet model. Note that the data and fit for B band were divided by a factor of 10, and that the data and fit for K' band were multiplied by a factor of 10 to 
avoid overlap between the three curves. 



